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FLUTTER OF BUCKLED, SIMPLY SUPPORTED, RECTANGULAR PANELS 

AT SUPERSONIC SPEEDS 

By Robert W. Fralich and John A. McElman 
Langley Research Center 

SUMMARY 

A theoretical flutter analysis is presented for buckled, simply supported panels 
subjected to supersonic flow over one surface. The analysis employs the Von Karman 
large-deflection plate theory and linearized static aerodynamic strip theory. A Galerkin 
procedure using four static mode shapes is employed to determine a set of differential 
equations which is programed on an analog computer. The character of the output of the 
analog is used to determine the flutter speed. Results are obtained for panels with ratios 
of length in the streamwise direction to length in the cross-flow direction equal to 1/2 
and 1 for three specified in-plane edge-loading conditions. An assessment of effects of 
cross- flow coupling of the modes is made by comparison of the results with those obtained 
when cross-flow coupling between the modes is neglected. 

INTRODUCTION 

Panel flutter has been encountered in the operation of aircraft and missiles and has 
become an important consideration in the design of structures for such vehicles. The 
panel flutter problem is influenced by many factors, such as aerodynamic effects, effects 
of boundary conditions and midplane compressive loads, and length-width ratio. In addi- 
tion, if the midplane compressive stresses are of sufficient magnitude to cause buckling, 
the flutter problem is further complicated. These various aspects of the flutter problem 
are discussed, for example, in references 1 to 10; some of the earlier investigations are 
listed in reference 1. Reference 1 considers the flutter of buckled, simply supported 
panels. In reference 1 the mechanism for flutter is shown to be a streamwise coupling 
between the modes. The purpose of the present analysis is to obtain a more accurate 
solution to this particular problem by investigating cross-flow coupling between the 
modes. Cross-flow coupling is a phenomenon which does not occur in the small- 
deflection flutter analysis of unbuckled, simply supported panels. 

The results presented in reference 1 were obtained from a two-mode Galerkin solu- 
tion, and rigorous analytical methods were used to determine the panel stability. In the 
present analysis a general Galerkin solution is obtained and the stability of the panel is 
determined for a four-mode solution by use of an analog computer. Points on the flutter 


boundaries have been determined for three in-plane edge-loading conditions for ratios 
of the length in the streamwise direction to the length in the cross-flow direction equal 
to 1/2 and 1. 

SYMBOLS 


AqjAijBq, 

*ljA m n»®mn 


n * D i coefficients in displacement expressions (eqs. (9)) 


Cmn 

D 

E 


length of plate in streamwise direction 


width of plate in cross-flow direction 


amplitude coefficients for lateral displacement 


flexural rigidity, 


Young's modulus 


Etr 


12(1 - p2) 


^mn>^mni^mn Fourier coefficients (see eqs. (12)) 


plate thickness 


M 


Mach number 


Mi, M2, M3, M4 

p l> p 2> p 3> p 4 

Q,S,T,J,K,N,H 


m,n 


parameters defined in equations (21) 
number of half-waves in streamwise and cross-flow directions, respectively 


NjjjNyjNjjy midplane stress resultants 

_ _ _ N x a 2 N y a 2 Nxya 2 

N x ,Ny,Nxy nondimensional midplane stress resultants; — - — , — - — , and J 


respectively 


7t 2 D ’ 7t 2 D : 


tt 2 D 


p x> p y average in-plane edge loads per unit length, positive in compression 


p xj p y 


nondimensional in-plane edge loads per unit length; 
respectively 


and 32-, 

7T 2 D tt2d 


2 


P,q,r,s 


integers 


q 

t 

u,v 

u,v 

V 

w 

w 

x,y,z 
£ = vm^ 

y 

®mn 

\ 

^•cr 

P 

y = a/b 
P 


1 9 

dynamic pressure, -pV* 

A 

time 

in-plane displacements, positive in x- and y-directions, respectively 

nondimensional in-plane displacements, positive in x- and y-directions 

u and v , respectively 

7T^D TT^Db 

free-stream velocity 

lateral deflection, positive in z-direction 
nondimensional lateral deflection, positive in z-direction, 
rectangular Cartesian coordinates (see fig. 1) 



mass density of plate material 


Kronecker delta, equals 1 if m = n, equals 0 if m ^ n 


speed parameter, 


16qa 3 


37t^/3D 

flutter speed parameter 


Poisson’s ratio 


nondimensional coordinates; x/a and y/a, respectively 


free-stream density of fluid 




nondimensional time, tl 

7 o 


2lW 


When subscripts £, i), and t follow a comma, they indicate partial differentia- 
tion with respect to £, 77 , and r, respectively. Dots over symbols denote derivatives 
with respect to r. 


STATEMENT OF PROBLEM 


y.n 


The configuration analyzed in this report is the simply supported, flat, rectangular 
panel shown in figure 1. The panel has a constant thickness h with air flowing over 

the top surface at a Mach number M. No flow of 
air beneath the panel is considered. Average 
in-plane edge loads P x and Py per unit length 
(positive in compression) are specified at the 
boundaries. No in-plane shearing forces are 
applied to the plate. 

METHOD OF SOLUTION 



Figure L- Rectangular panel and coordinate 
system. 


The present analysis employs the large- 
x, 5 deflection plate theory of Von Karman and linear- 
ized static aerodynamic strip theory. This aero- 
dynamic approximation has previously been shown 
to yield accurate flutter boundaries for Mach num- 
bers greater than about 1.6. The resulting equa- 
tions are analyzed in the appendix by means of a 
Galerkin procedure which utilizes the doubly 
infinite set of static buckling modes. This pro- 
cedure yields a doubly infinite set of second- order nonlinear ordinary differential equa- 
tions for the time- dependent amplitude coefficients C mn . These equations can be 
reduced to those for various approximate analyses that utilize a finite number of static 
buckling modes. The set of four equations, obtained from an approximate analysis that 
uses four of the static modes, is analyzed by means of an analog computer in order to 
find a flutter speed parameter Xcr- The modes considered have amplitude coefficients 
C mn , where m = 1 and 2 and n = 1 and 2 are the number of half waves in the stream - 
wise and cross-flow directions, respectively. 


In the analog analysis, the character of the time histories of the amplitude coeffi- 
cients is observed while the initial static buckling condition is altered by gradually 
increasing the speed parameter X to a given level. Analog traces which illustrate the 
method of determining flutter are given in figure 2. For levels of X below the critical 
value X cr (fig. 2(a)), the amplitude coefficients do not build up with time and the motion 
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(a) Stable motion. (\ < 
Figure 2.- 


A cr ) (b) Unstable motion. (A > A cr ) 

Variation of amplitude coefficients with time. 


is considered stable. For levels of X above the critical flutter speed X cr (fig. 2(b)), 
the amplitude coefficients show a drastic buildup with time. The character of the motion 
is observed for different levels of X until the critical flutter speed parameter X cr is 
determined. The value of X cr for the case illustrated in figures 2(a) and 2(b) is indi- 
cated by the tick mark. Values of X cr are determined for three in-plane loading con- 
ditions for two values of a/b, the ratio of length in the streamwise direction to width in 
the cross-flow direction. 


RESULTS AND DISCUSSION 

Results for the flutter boundaries obtained from the present analysis are given by 
the circles in figures 3 to 5. Flutter boundaries (x = X C r) are plotted in these figures 
as a function of in-plane edge loading for a/b of 1/2 and 1. Flutter occurs above 
these boundaries and stable motion is obtained below the boundaries. Results for speci- 
fied values of P x with Py = 0, for specified Py with P x = 0, and for specified 
P x = Py are shown in figures 3, 4, and 5, respectively. Also shown in these figures 
are the flutter boundaries presented in reference 1 which were obtained by a rigorous 
stability analysis using just two modes. The curve labeled n = 1 is the flutter 
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boundary obtained by considering only the Cu 

2 modes (unbuckled) ) , f , > and Coi modes, whereas the one labeled n = 2 

2 modes (buckled) )' re ' ' .. 

o 4 modes (present analysis) utilizes the C^2 an d C22 modes. 



(a) v = l/2. 



0 5 10 15 20 25 30 


(b) v = 1. 

Figure 3.- Flutter boundaries for pan_els with 
streamwise compressive load P x . Py = 0. 


The results show effects of cross-flow 
coupling of the modes on the flutter of simply 
supported, rectangular panels in supersonic flow. 
Theoretically, this coupling does not exist on the 
portion of the flutter boundary prior to buckling 
(the dashed line portions of figs. 3, 4, and 5). One 
result in the present four- mode analysis was 
determined in the unbuckled range. (See fig. 4.) 
Figure 4 confirms the absence of cross-flow 
coupling in the unbuckled range and gives an 

2 modes (unbuckled) ) < ref> ± ) 

2 modes (buckled) f 

o 4 modes ( present analysis ) 



2 modes (unbuckled) 1 , f n 

2 modes (buckled) f 

o 4 modes ( present analysis ) 



(a) v = 1/2. 



(b) v = 1. 


Figure 4.- Flutter boundaries for panels with 
cross-flow compressive load P y . P x = 0; 
v = 1. 


Figure 5.- Flutter boundaries for panels 
with equal streamwise and cross-flow 
compressive loads. P x = Py. 
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indication of the validity of the procedure used in determining the flutter criterion. On 
the portion of the flutter boundary for buckled panels, cross-flow coupling of the modes 
has an effect on the flutter boundary for certain values of a/b and in-plane loading 
conditions. For the two square panels under streamwise loading and under biaxial 
loading (figs. 3(b) and 5(b)) the effects are negligible. However, for the panels shown in 
figures 3(a), 4, and 5(a), effects of cross-flow coupling on the flutter boundary are sig- 
nificant. For these cases (except for a region of fig. 3(a)), the effect of cross-flow 
coupling is to lower the flutter boundaries. 

The results show that the effects of cross-flow coupling can be important in deter- 
mining flutter boundaries for buckled panels. In this analysis no investigation is made 
of the convergence of the Galerkin solution. To do this more modes would have to be 
considered. In order to facilitate such a study, the modal solution that includes all the 
static buckling modes is presented in the appendix. 

CONCLUDING REMARKS 

A supersonic flutter analysis is presented for simply supported, rectangular panels 
subjected to specified in-plane compressive edge loads. A Galerkin procedure that 
utilizes the doubly infinite set of the static buckling modes yields a doubly infinite set of 
differential equations that can be reduced to any desired finite number of equations. The 
equations have been programed on an analog computer for a stability analysis for a four- 
mode solution that exhibits both streamwise and cross-flow coupling of the modes. The 
character of the output of the analog computer is then used to determine the critical 
flutter condition. Numerical results are presented for panels with ratios of length in the 
streamwise direction to length in the cross-flow direction equal to 1/2 and 1. The fol- 
lowing specified in-plane boundary edge conditions are considered: (a) streamwise com- 

pressive loading only, (b) cross-flow compressive loading only, and (c) equal streamwise 
and cross-flow compressive loading. The results show that the effects of cross-flow 
coupling are important in determining flutter boundaries for buckled panels since the 
flutter speed can be appreciably reduced from the value determined without cross-flow 
coupling. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 6, 1967, 

126-14-02-24-23. 
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APPENDIX 


ANALYSIS 


The nonlinear differential equations expressing the equilibrium of an aerodynami- 
cally loaded oscillating panel based on Von Karman large-deflection plate theory are 
obtained in nondimensional form in reference 1 . These equations can be written as 
follows: 


( 1 ) 


+ i y ii{ \w + i f J£ a ^2 ( W , + MW ,1? 2 ),4 " ^2 

a + + = ■ + ^^ 2 ), n '^2 (2) 


+ 2w ,^66 + W ,* 7 W?) - ^( N x w ,^ + + 2N xy W , $7) + I Xw ,£ + w ,rr = 0 


3 


(3) 


where 


1 

-^ 2 

L- 


1 

[b 

1 

-M 2 

L a 


1 

1 


Mv 2 -= 2 


,« +mw j?7 


b - - 1 /— 2 - 2\ 




(4) 


in which u, v, and w are nondimensional displacements, N x , Ny, and Nxy are non- 
dimensional stress resultants, £ and rj are nondimensional x and y coordinates, 
r is nondimensional time, j u is Poisson's ratio, and A is the speed parameter. 

The boundary conditions to be satisfied by w are the simple-support conditions 

-N 

w( 0,?7 ,t) = w(l ,?7 ,r) = w(^,0,t) = = 0 

w ^( 0 ,v,r) =W f ^(l,?7,T) =w >T77 j(|,0 ,t) = w >7 ?T7 (^,t) =0 




(5) 


The boundary conditions to be satisfied by u and v are those for uniform displace- 
ment of each edge in the plane of the plate 

u 77(0,77,7) =u 77(1,77,7) = v ^(£,0, t) =v^,jt) = 0 


( 6 ) 
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appendix 


and those for zero in-plane shear stress at the edge of the plate 

v |(0,?7 ,t) = v^(1,77,t) = u >7? (£,0,t) ^(^t) = 0 

Appropriate boundary conditions for the edge loads are 

f I 0 b/a (^){=0 - 5 1 0 b / a ( H = i d ” - 

J 0 ( N y)r?=o = j * 0 ( %) 77= b d| = - p y 


( 7 ) 




( 8 ) 


and P x and Py are nondimensional in-plane edge loads per unit length. 

The boundary conditions (eqs. (5), (6), and (7)) can be satisfied if the displacements 
u and v are written as 


( 9 ) 



OO OO 


u = Aq + A^ + 

I I 

A mn sin m7r ^ cos 


m=l n=0 


v = Bq + B-^ + 

00 oo 

1 l 

B mn cos mir £ sin 


m=0 n=l 

J 


and if the normal displacement w is written as 

OO OO 

w 

m=l n=l 

The nonlinear terms on the right-hand side of equations (1) and (2) can also be 
expanded in Fourier series as follows: 


oo oo 

- ^ ^ Cmn sin sin n7T ^ r l 


( 10 ) 


1 

27 



^ ^ F mn cos cos 

m=0 n=0 



£ £ G mn cos cos MZ > 

m=0 n=0 


OO OO 


3 w ,« w ,» = I I H mn sin sin SSS 

It 


m=l n=l 


(id 
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where 


^mn _ o 

b7T“ 


f f w t 2 COS mir^ cos 2~Z d£ d?7 

( 1+5 mo)( 1 + 6 0") J ° J ° 

f f 5 _ 2 cos mm? cos MM d| djj > 

( 1 + «mo)( 1 + 5 no) J ° J » b I 


'mn 9 
bir^ 


(12) 


Hmn - ^5 c i; w ,5 w ,n sin m ’ s sin tt 2 d « 


drj 


where 


6 mn ~ 0 
®mn = * 


(m * n) 
(m = n) 


When equations (9) and (11) are substituted into equations (1) and (2) the coefficients 
A mn and B mn can be determined in terms of F mn , G mn , and H mn as follows: 


1 mn - 




- -\m 
2 ir 


m 2 + 


<*♦“>(?) 


F mn + 


jum‘ 


-(f) “ G — (?) H 


mn 


(m = 1, 2, 3, n = 0, 1, 2, . . ■*>) 


Bmn = ~ 


m^ + 


(?) 


r 

r- 



” 

-j 



na 

l b 

(2 + /z) m 2 + 0 

ft 

G mn + 

if 

f- m 2 

h 

11 

pj F mn + mH mn j 


(m = 0, 1, 2, . . <*>; n = 1, 2, 3, . . «>) 


J 


(13) 


The displacements u and v are now known in terms of Ag, Aj_, Bq, Bj_, 
F mn , G mn , and H mn . The stress resultants are obtained in terms of the coeffi- 
cients Aj, Bj, F mn , G mn , and H mn by substitution from equations (9), (11), and 
(13) into equation (4). In order for N x and Ny to satisfy the boundary conditions 
(eq. (8)), the constants Aj and Bj must have the values 


A 1 - _p x + M p y ‘ F 00 
®1 = g(- p y + ^ p x " Gggj 




(14) 


j 
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where, when use is made of equations (10) and (12), 


oo CO 


F o° = l Z I 

p=l q=l 


G °° “ I ^2 Z Jj ( l 2c pq 2 

b p=l q=l 




( 15 ) 


The stress resultants become 


N X = -P X+ ; F 0n cos22S + v 

m= 


I 

n=l 


(ft 


J mn 


- - [ m * ♦ (ff 


cos mir£, cos 


b 


Ny = -Py + 


oo oo oo 

^ G m0 cos miri; + I I 


m 2 J, 


mn 


s *y = I Z 


m=l 
~ m 


m=1 n=1 m 2 + f — V 

\ b / 


cos mni; cos MHZ 
* b 


> (16) 


mn 


m=l n=l 


m 2 + 


{f)‘ 


sin m7r| sin n — 1 


where 


Jmn ~ ("b - ) ^mn + m ^mn + (11) 

When use is made of equations (12) and (14), F 0n , G m Q, and J mn become 


oo oo oo 


F 0n - | X X Jy p2c Pq C Ps[ 26 (s-q)n ' 6 (q+s)^] 
p=l q=l s=l L - 1 


OO 00 OO 


G m° “ f p 2 Jy /L q2c pq C rq[ 26 (r-p)m ~ 6 (p+r)m] 


p-1 r=l q=l 


(n = 1, 2, . . °°) 


(m = 1, 2, . . °°) 


} (18) 


g OO OO OO 03 ~ 

Jmn = I ^2 S I 1 1 C pq C rs{ 2 ( mc l + np)(ms + nr)6 (r _ p)m 6 (q _ s)n + (mq - np)(ms - nr)[26 (r . p)m 6 (s . q)n 
D p=l r=l q=l s=l v 

* Vr)m 6 (q + s)n] + 2(mq + nP)(mS ' nr) [ 6 (r-p)m 6 (q + s)n “ 6 (p +r )m 6 (s-q)n]\ 
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Substitution of equations (10) and (16) into equation (3) and application of the Galerkin 
procedure yields 


d C mn 
dr^ 




m 2 P x 





2 00 OO er\ 

^ 2 °° °° 

+ TZ Z C msF 0q [«( s+q )n + 6 (s-q)n ' 6 (q-s)n] + £(^) £ £ C m G p0p( r+ p)m + 6 (r-p)m " 5 (p-r)m] 

S_1 q ~ 1 r=l p=l 


_ a 2 V V V V C rsJpq r 

+ t-tt ) ) ) ) ~ =j2 |_ 6 (r+p)m 5 (s+q)n + 6 (r-p)m 6 (s-q)n ” 5 (r-p)iri (q-s)n 


■mm- 

D r=l p=l s=l q=l 


P 2 + 


/gar 

Vb ) 


OO OO OO OO 


i 1 a 2 V V V V ^ rq + SP ^ Cr s J pq r , 
" 6 (p-r)m 6 (s-ci)n + 6 (p-r)m 5 ( q -s)nJ + 4 ^2 L L L L, ~ F " ~f|2 \_ 

J b r=l p=l s=l q =l p 2 + m 2 

6 (r+p)m 5 ( q -s)n + 6 (r-p)m 6 (s+ q )n 5 (p-r)m^(s+ q )n] “ 0 


6 (r+p)m 6 (s- q )n 


( 19 ) 


The modal solution (eq. (19)) that utilizes the doubly infinite set of static mode 
shapes reduces to any desired approximate solution that uses a finite number of modes 
simply by deleting the undesired modes. In the present analysis numerical results are 
found for a four-mode solution so that 

w = sin + C 21 sin 27r^sin sin + C 22 sin 27r^jsin 2^2 (2 0) 

and equation (19) becomes 

Cn + c n (- Mj + PiC n 2 + QC 2 1 2 + SC 12 2 + TC22 2 ) + C2l(HC 12 C 2 2 " a) = 0 

C 21 + c 2l(“ m 2 + Q°ll 2 + P 2 C 21 2 + JC 12 2 + KC 22 2 ) + Cll( HC 12 G 22 + a) = 0 
c 12 + c 12 (- M 3 + sc n 2 + JC 21 2 + P 3 C 12 2 + NC 22 2 ) + C 22 (HC n C 21 - x) = 0 
^22 + ^ 22 ^ ^4 + TC-ll 2 + KC 2 j 2 + NC 12 2 + P4.C 22 j + Ci2(HCnC2i + a) = 0 


where 


2 

Ml = P X + P 2 Py - (l + 1/2) 

Q 

M 2 = 4P X + J/ 2 Py ~ (4 + 1/2) 


( 21 ) 
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M 3 = P x + 4y2py - (l + 

M 4 = 4P X + 4^ 2 Py - 16(1 + y2 ) 2 

p i-'ti( 14 ‘- V4 ) 

p 2 = + p4) 

P 3 -^(l+ 16 ^ 4 ) 

P 4 = 1 + 


Q = 


16 


4(1 + xA) 


8I1A x 1 A 

(l + 4 (9 + 4,2) 2 


s = — 

16 


T = 16 


4(l + v 4 ) 


81f/ 4 , i/ 4 

H H 


(4 + z/2) 2 (4 + 9y 2 ) 


(l + 9z/ 2 ) (9 + z/ 2 ) 


1 

62 5i/ 4 

16 

(1 + 9,2 f 


t 625fr> 4 | 82^ 4 


K = 


N = 


jtw,/*)— SIS* ♦ 


(l6 + V 2 ) 2 (16 + °>v 2 f 


i(l + 16*4) + 81y l - + £— 

(l + 16 v 2 ) 2 (9 + 16^ 


H = 1 + i/ 4 + 25 ^ 


25^ 


(l + 9^) 2 ( 9 + v 2 ) 


in which v = a/b. In equation (21) the double dots over the coefficients C m n repre- 
sent the second derivative with respect to r. Equations (21) are set up on an analog com- 
puter in order to determine the flutter speed X = \ cr as a function of the in-plane edge 
loads and a/b. 
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